A Geometrical Approach to Strong Gravitational Lensing in f(R) Gravity 
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We present a framework for the study of lensing in spherically symmetric spacetimes within 
the context of f(R) gravity. Equations for the propagation of null geodesies, together with an 
expression for the bending angle are derived for any f(R) theory and then applied to an exact 
spherically symmetric solution of R n gravity. We find that for this case more bending is expected 
for R n gravity theories in comparison to GR and is dependent on the value of n and the value of 
distance of closest approach of the incident null geodesic. 



I. INTRODUCTION 

Modifications of the theory of General Relativity (GR) 
by including higher order corrections to the theory were 
first introduced by Herman Weyl [l[ as early as 1918, 
just three years after Einstein completed his theory of 
General Relativity. Weyl proposed modifications to the 
theory by including higher order invariants in the action, 
in order to extend general relativity's geometrical founda- 
tions. In 1921, Arthur Eddington also began to consider 
fourth order theories of gravity and went on to publish 
his famous book: The Mathematical Theory of Relativ- 
ity which contained work on generalized versions of 
Weyl's theory. Ever since, there has been a great number 
of gravitational theories proposing modifications to the 
Einstcin-Hilbert action. 

Originally the motivation for considering modified the- 
ories of gravity stems from limitations of General Rel- 
ativity when considering strong gravity regimes. More 
recently however, corrections to GR have been proposed 
to account for the dark sector of the universe. Studies 
of the Cosmic Microwave Background (CMB) Q, Su- 
pernovae Type la surveys [ij and Baryon Acoustic Os- 
cillations |5j indicate that the energy density budget of 
the universe is 5% ordinary matter (baryons, radiation 
and neutrinos), 25% dark matter and 70% dark energy. 
Dark matter is responsible for the gravitational clump- 
ing of galaxies, galaxy clusters and large-scale structures, 
whilst dark energy is linked to the current phase of ac- 
celerated expansion of the Universe. The ACDM model 
(concordance model) is currently the best fit model, us- 
ing the cosmological constant A to explain dark energy. 
An alternative explanation is to consider that on cosmo- 
logical scales, the current description of gravity is incom- 
plete and a modification to the theory is needed. One 
of the most popular modifications to gravity which has 
been proposed to explain the late-time acceleration of the 
Universe, without the need for dark energy, is f(R) grav- 
ity, which represents a generalization of Einstein-Hilbert 
action by replacing the Ricci scalar R with the function 



f{R). These fourth order theories have been shown to 
and give rise to a phase of accelerated expansion with- 
out the need for dark energy |6|-|9j and to account for the 
rotation curves for spiral galaxies without the need for 
dark matter [l(| ■ 

A good theory of gravity is clearly one that agrees 
with observations. One of the first experimental veri- 
fications of GR was the measurement of the bending of 
light, observed for the first time by Eddington in 1919 
during a total solar eclipse. Since then, Gravitational 
lensing has been a powerful tool used to determine the 
mass distribution of galaxies and galaxy clusters and to 
put constraints on scales as small as stars to large scale 
structures and cosmological parameters [ll[. Given that 
the lensing effect is dependent on the underlying theory 
of gravity, investigating modifications of GR would re- 
sult in deviations from the standard expression of the 
deflection angle and is consequently worth investigating. 
The derivation of the form of the lensing angle for f(R) 
theory has already been presented in [l2j and [l3| using 
different methods. 

In this paper we use covariant method to find general 
expressions for the propagation of null geodesies in static 
spherically symmetric spacetimes in f(R) gravity theo- 
ries. The solutions of these equations are then used to 
obtain a general expression for the deflection angle for 
spherically symmetric spacetimes. The outline of this 
paper is as follows. In section [TT] we present the gen- 
eral equations for f(R) gravity. In section HTT1 we apply 
thel+1+2 covariant approach to static spherical sym- 
metric spacetimes in f(R) gravity. The propagation of 
null geodesies is considered in section IIVI and then cal- 
culated explicitly for the case of f(R) = R n in static 
spherically symmetric spacetime in order to determine 
the bending angle in section IVl 

Unless otherwise specified, natural units (h = c = 
ks = 8ttG = 1) will be used throughout this paper, Latin 
indices run from to 3. The symbol V represents the 
usual covariant derivative and d corresponds to partial 
differentiation. We use the —,+,+,+ signature and the 
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Riemann tensor is denned by 

n bed = 1 bd,c — i bc,d + -L bd-L ce ~ t bet de , {*-) 

where the T a bd are the Christoffel symbols (i.e. symmet- 
ric in the lower indices), defined by 

T Q M = ^g ae {g be ,d + 9ed,b - g bd} e) ■ (2) 

The Ricci tensor is obtained by contracting the first and 
the third indices 

Rab — g Cd Racbd ■ (3) 

The symmetrization and the antisymmetrization over the 
indexes of a tensor are defined as 

T( a b) = 2 ( Tab +Tba) , T[ ab ] = - (T a b -Tba) ■ (4) 

Finally the Hilbert-Einstein action in the presence of 
matter is given by 

A =\f <£xj—g [R + 2C m ] . (5) 



II. GENERAL EQUATIONS FOR FOURTH 
ORDER GRAVITY 

In a completely general context, a fourth order the- 
ory of gravity is obtained by adding the terms involv- 
ing f{R), RabR ab , RabcdR abcd to the standard Einstein 
Hilbert action. However, we know the Gauss-Bonnet 
term (Q = R 2 - 4R ab R ab + R a bcdR abcd ) is a total dif- 
ferential in four dimensions and hence do not affect the 
field equations. Using this symmetry we can rewrite 
R a bcdR abcd in terms of the other two. Furthermore, if 
the spacetime is highly symmetric, then the variation of 
the term R a bR ab can always be rewritten in terms of the 
variation of R 2 [3, [ill- It follows that the sufficiently 
general fourth-order Lagrangian for a highly symmetric 
spacetime only contains powers of R and we can, with 
out loss of generality, write the action as 

A=\J<fx^—g[f(R) + C m ], (6) 

where C m represents the matter contribution and A is 
the usual cosmological constant. Varying the action with 
respect to the metric gives the following field equations: 

f,R G ab =T^ + \{f- Rf,R ) + V b V a f,R - 5 afcV c V C /,fl , 

where denotes the derivative of the function '/' w.r.t 
the Ricci scalar and T"? is the matter stress energy tensor 
defined as 

T™ = fi m U a U b + p m h ab + Cu b + qTu a + TT™ b . (8) 



Here u a is the direction of a timelike observer, h a b is 
the projected metric on the 3-space perpendicular to u a . 
Also /j, m , p m , q m and 7r™£ denotes the standard matter 
density, pressure, heat flux and anisotropic stress respec- 
tively. Equations (J7]) reduce to the standard Einstein 
field equations when f(R) = R. 



III. 1+1+2 COVARIANT APPROACH 
A. Relativistic Cosmology 

We know that the 1+3 covariant approach, initially de- 
veloped by Ehlers and Ellis l (j has proven to be a very 
useful technique in many aspects of relativistic cosmol- 
ogy, In cosmology these methods have been applied to 
the formalism and evolution of density perturbations in 
the universe 1171 and to the physics of cosmic microwave 
background [18j amongst other things. This approach is 
based on a 1+3 threading decomposition of the space- 
time manifold w.r.t a timelike congruence as a splitting 
of spacetime onto a timelike and an orthogonal three- 
dimensional spacelike hypersurface. All the essential in- 
formation in the system is captured in a set of kinematic 
and dynamic 1+3 variables that have a well defined phys- 
ical and geometrical significance and that satisfy a set 
of evolution and constraint equations derived from the 
Bianchi and Ricci identities, forming a closed system of 
equations for a chosen equation of state describing mat- 
ter. 

A natural extension to the 1+3 approach, optimized for 
problems which have spherical symmetry, is the 1+1+2 
formalism developed recently by Clarkson and Barrett 
[l!| . In this formalism, a further splitting of the spacelike 
1+3 variables is performed such that it isolates a specific 
spatial direction. This allows for the derivation of a set 
of variables that are more advantageous to treat systems 
with one preferred direction. For example in spherically 
symmetric system the equation for the 1+1+2 variables 
are scalar equations and are much simpler than the ones 
of the 1+3 formalism which are in general tensorial. The 
1+1+2 formalism was applied to the study of linear per- 
turbations of a Schwarzschild spacetime [19| and to the 
generation of electromagnetic radiation by gravitational 
waves interacting with a strong magnetic field around a 
vibrating Schwarzschild black hole |20| . 

In (1+3) approach first we define a timelike congruence 
by a timelike unit vector u a . Then the spacetime is split 
in the form R ® V where R denotes the timeline along u a 
and V is the 3-space perpendicular to u a . Then obviously 
u a u a = — 1, and any vector X a can be projected on the 
3-space by the projection tensor h b = g b + u a Ub- 

Two derivatives are defined: the vector u a is used to 
define the covariant time derivative (denoted by a dot) 
for any tensor T c ^ a " b along the observers' worldlines de- 
fined by 

T c .. d a -» = u e V e T c .. d a -» , (9) 
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and the tensor h a b is used to define the fully orthogonally 
projected covariant derivative D for any tensor T c ..d a ' , 

D e T c .. d a - b = h c P..h d ih f a ..h g b h/V r T p .. q f-9 , (10) 

with total projection on all the free indices. In the 
(1+1+2) we further split the 3-space V of the (1+3) ap- 
proach, by introducing the unit vector e a orthogonal to 
u a so that 



e a u a = , 



= 1. 



(11) 



effective "fluid" arising form a modification of the un- 
derlying theory of gravity. We will see that this will be 
actually the case in the context of fourth order gravity. 

The full pre sentation of the 1+1+2 equations can be 
seen in Halm. 



B. Spherically symmetric static spacetimes in f(R) 
Gravity 



Then the projection tensor 

N a b = h a b - e a e b = g a b + u a u b - e a e b , N\ = 2 , (12) 

projects vectors onto the 2-surfaces orthogonal to e a and 
u a , which we would refer as 'sheets'. 

Apart from the 'time 1 (dot) derivative, of an object 
(scalar, vector or tensor) which is the derivative along 
the timelike congruence u a , we now introduce two new 
derivatives, which e a defines, for any object ip a ...b c '" d - 



i>a. 



c.d 



c.d 



e f Dfip a .. b 



c.d 



(13) 



Nj ...N b 9 N h c ..Ni d Nf D^fJ~i .(14) 



The hat-derivative is the derivative along the e a vector- 
field in the surfaces orthogonal to u a . The ID-derivative is 
the projected derivative onto the sheet, with the projec- 
tion on every free index. The irreducible set of kinematic 
and dynamical quantities: 

[e, <j>, Cab, A, A a , a a , a a , & n, n a , (15) 

E, E Q , E ab , £, £ a , £ ab , H, U a , U ab ] (16) 

are obtained from the 1+1+2 splitting of the Weyl ten- 
sors, the covariant variation of e a and the Ricci and 
Bianchi identities for e a and u a . 

We see that travelling along u a , is the expansion 
scalar (volume expansion), whilst travelling along e a , 4> 
represents expansion of the sheet and Cab is the shear of 
e° (i.e., the distortion of the sheet). A is the radial com- 
ponent of the acceleration of u a , that is, of ii a [27j and 
A a is its component lying in the sheet orthogonal to e a . 
a a is the component of n a lying in the sheet and a a is 
the acceleration of e a . £ can be interpreted as the vortic- 
ity associated with e a , il and il a are the components of 
the 3-vector vorticity uj a along e a and its sheet compo- 
nent respectively. The scalar, vector and tensor parts of 
the projected, symmetric, trace-free of the 3-tensors are 
E , E a and E ao , respectively, for the shear, £ ,£ a ,£ a b for 
the electric Weyl tensor and H , H a , H a b for the magnetic 
Weyl tensor. 

The quantities (IT5l) . together with the set of thermo- 
dynamic variables 



[^p,Q,U,Q a ,U ab 



(17) 



associated to a generic fluid present in the system, com- 
pletely characterize the 1+1+2 spacetime. These ther- 
modynamic quantities can be also used to represent the 



LRS spacetimes posses continuous isotropy group at 
each point and hence a multi-transitive isometry group 
acting on the spacetime manifold 22] . These spacetimes 
exhibit locally (at each point) a unique preferred spatial 
direction, covariantly defined, for example, by either vor- 
ticity vector field or a non- vanishing non-gravitational ac- 
celeration of the matter fluids. The 1+1+2 formalism is 
therefore ideally suited for covariant description of these 
spacetimes, yielding a complete derivation in terms of in- 
variant scalar quantities that have physical or direct geo- 
metrical meaning (23| . The preferred spatial direction in 
the LRS spacetimes constitutes a local axis of symmetry 
and in this case e a is just a vector pointing along the axis 
of symmetry and is thus called a 'radial' vector. Since 
LRS spacetimes are constructed to be locally isotropic, 
this allows for the vanishing of all 1+1+2 vectors and 
tensors, such that there are no preferred directions in 
the sheet. Thus, all the non-zero 1+1+2 variables are 
covariantly defined scalars. The variables, 



H, e, & £,E,n, p,n,Q] 



(18) 



fully describe LRS spacetimes and are what is solved for 
in the 1+1+2 approach. A detailed discussion of the 
covariant approach to LRS perfect fluid spacetimes can 
be found in (Hill. 

The LRS class II spacetimes admit spherically sym- 
metric solutions that are rotation free. From this and 
the propagation and commutation equations characteriz- 
ing this space time, the variables f2, £ and % vanish [23| . 
The condition of staticity further implies that the dot 
derivatives of all the quantities vanish and consequently 
E, G and Q also vanish. 

These conditions reduce the full LRS system of equa- 
tions to: 



^2 2 



1 



-u--n-£ 

3 P 2 







P + n 
A 



- 3 A £+ l n 



(19) 
(20) 



Acf> + ^{ l i + 2 > p)-£ + in ,(21) 

3 
2* 



A)Il-( f i + p)A, (22) 
1 



(A + ^A+^ifi + Zp) 



(23) 
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The quantities fi, p and II are denned, in this 



M = ~J~ (V" + l(Rf,R -/) + f,RR X 
J,R V 1 



~f,RR X(f) + f,RRR X 2 ) 



P = -J- (p m + kf~ Rf,R ) - h,RR X- 

j,r \ z •j 
2 2 

-^f,RRX(/)—-f,RRRX -Af,RRX 



x a (v), where v is an affine parameter along the geodesies. 
The components of the tangent vector of x a (v) are 



(24) 
(25) 



dr a 

and it obeys: 

fc a fc Q = , 
From this, it follows that 



n 



1 /2 2 2 "1 

\^f,RRR X + ~f,RR X - -f,RR X(f> 



(26) 



(27) 



Using the equations in (|19H23p including (|28|) and elim- 
inating £, we get the set of four coupled first order equa- 
tions governing the spacetime in the fourth order gravity 
as 



f,R 



f,R 



A + A{A - 



\ Rf,R -\ f 
+f, RR X{cf> + 2A) , 



(28) 



1 



1 



R = 



J,rrX — 



G 
X 



f - TT^/^i? fjRR XA , 



(29) 



(30) 



1 

-XI 



—■^Rf,R 



f - f,RRR X z 



3 

A)f,RR 



(31) 



From the system of equations (|28H31j) we can deduce 
some important results for spherically symmetric static 
solutions in a general f(R) gravity in a completely co- 
ordinate independent manner (see [24j]). 



k b V b k a 



Sk a 



8v 
s _ /> 



= , 



(32) 



(33) 



(34) 



where we have defined 4- — k Vb as the derivative along 
the ray. Since k a is a gradient, we have that Vbk a = W a k b 
and this implies that (|34j) can be written as 



where we have defined R = X . We will consider the de- 
flection angle of light propagating near an isolated spher- 
ically symmetric mass distribution so that /x m = and 
p m = 0. Because of the additional degrees of freedom the 
equations (|19H23[) are not closed and we have to add the 
trace equation to achieve closure: 

RflR—'Zf = -3f,RRX — 3f,RRX(f) — 3f,RRRX 2 — 3Af,RRX 



k b W a k b = . 



(35) 



The light propagation vector k a is received by the ob- 
server from the direction determined by the the unit spa- 
tial vector n a : 



1 







(36) 



The null vector k a can be split into in the usual way [ll| : 
k a = E(u a + n a ) , (37) 
A 1 + 1 + 2 split of n a can then be performed (2f| , giving 



k a = E{u a 



ne 



(38) 



where E = —u a k a and since u a k a is the circular fre- 
quency of the electromagnetic wave, E can be interpreted 
as the energy associated with the ray. k is the magnitude 
of the radial component and K a is the component lying in 
the 2-dimensional sheet. Using this form of the null vec- 
tor, we can determine the lensing geometry of a photon 
experiencing a deflection about the centre of symmetry. 

In the presence of a strong gravitational field such as 
a black hole, with k a lying tangent to the null geodesic, 
the general scalar deflection angle takes the form [251 ] : 



f 2 - \E\ ^J\-K 2 dv- 



a 



(39) 



This relation is geometrical and completely general. As 
a consequence, if we know the quantities E(u), n{v) and 
r(u), for a given spherically symmetric spacetime, it is 
possible to find an explicit form of the deflection angle. 



IV. LENSING GEOMETRY 

A. Null geodesies 

We now apply the 1+1+2 approach to null geodesies 
characterized by a the family of null curves (or light rays), 



The propagation equations for the lensing 
variables. 



Let us now look at the general propagation equations 
for the lensing variables E and k in the direction of the 
null ray k a . The geodesic condition ([3"5"j) can be used to 
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derive propagation equations for E and k. Substituting 
for the null vector k a into (|55|) and projecting the 
expression along the timelike direction (u a ) and along the 
radial direction (e a ), we obtain the general propagation 
equations for E and k: 

§ = ^ = - E ^-| s ^-^(ie-i E ),(4o) 



coordinate r. The most natural way to do this is to make 
the Gaussian curvature 'K' of the spherical sheets to be 
proportional to the inverse square of the radius K = 
■K. In that case, this co-ordinate V becomes the area 
radius of the sheets. This gives a geometrical definition 
to the 'hat? derivative. From the propagation equation 
K, K — —<j>K [23l ]. the most natural way to define the 
hat derivative of any scalar "J/ would be 



8 k 



E (l 



(41) 



respectively. Here, spherical symmetry has been consid- 
ered and the properties 

k b u b = -E, k b e b = En, N a b k b = En a , 
e a b k b = Ee a b K b , u a n a = 0, e a n a = (42) 

have been used. We have also utilized the expressions 



u' a = EAe a + EK + S 1 ' 

+E Qe- J K a + EVte ab n b 



(43) 



EAa„ Eh j S + -6 j //„ 



1 



-S(/>K a + E£e ab n b , 



which are obtained from (31) and (33). In LRS space- 
times, the general propagation equations (|4"01) and (HTl) . 
in the direction of the ray, reduces to 



E' = -E 2 Ak , 

k' = E(1-k 2 )(^-A) . 



(44) 
(45) 



It is worth noting at this point that the form of the above 
equations is independent of the theory of gravity. As con- 
sequence differences in the lensing phenomenon between 
GR and other theories of gravity will only depend on the 
features of the metric of the gravitational source. 

When considering spherical symmetry, it is convenient 
to express the hat derivative in terms of the proper radial 



t 1 



(46) 



for a static case. Using this we can write the propaga- 
tion equations for E (|4"4")) and k (|43|) in terms of the hat 
derivative as 

E' = k a V a E = EkE , (47) 
k' = k a V a K = Ekk. (48) 
Expressing the hat-derivative in terms of the radial pa- 
rameter r and rearranging, we obtain: 



(1 



1 dE 
E~dr 

K dhl 

— k 2 ) dr 



2A 

r<f> 



1 



2A 

r<p 



(49) 
(50) 



In order to find solutions to these differential equations, 
one has to first obtain solutions to the system (I28H31[) to 
obtain general forms for A,<p ,R and X and then substi- 
tute them into (gSJ) and (|50]l. Once the solutions of (gSJ) 
and are found, one can substitute them into equa- 
tion (|39|) to calculate the deflection angle in the presence 
of a strong gravitational field. 



V. A SIMPLE EXAMPLE: iT-GRAVITY 

A. Solution for the lensing variables 

Let us now consider the simple case in which f(R) = 
\R n • We will use a solution for the quantities A,4>,R 
and X found in [U]. They read: 
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A = 



-C (5 - An) r 4 " "-^" +9 + (4n 2 - 6n + 2)r~ 1 
2(2-ra) 



(l + 2n - 2n 2 ) (7 - lOn + 4n 2 ) ( 1 + Cr * 
(2-n) 2 



fl + 2n - 2n 2 ) (7 - lOn + 4n 2 ) 



V 



i? = 



X = 



(2-n) 2 ( 
6n(n — 1) 



1 + CV 



(2n (n-1)- l)r 2 ' 

12n (n- 1) | (1 + 2n - 2n 2 ) (7 - lOn + 4n 2 ) 
"(2n(«-l)-l)r3 ^ (2 _ n)2 ^ 1 + Cr -l-^i^ 



(51) 

(52) 
(53) 

(54) 



It is important to remember that this metric can be used 
only for n < (1 + y/o)/2 m 1.23. Beyond this value of 
n the signatures changes and the solution above must 
should be considered unphysical in this context. It is 
worth noting here the generalization of Birkhoff 's Theo- 
rem obtained in [24j for higher order gravity states that 
for all functions f(R) which are of class C 3 at R = and 
/(0) = while f,R (0) ^ 0, the Schwarzschild solution 
is the only static spherically symmetric vacuum solution 
with vanishing Ricci scalar. This effectively implies that 
the form of the deflection angle for this case would be 
the form of the angle obtained in (2(| for a Schwarzschild 
spacetime and the standard lensing results will be ob- 
tained. 

Now substituting for A and 4> from (l5Tj) and (|52p re- 
spectively in (|49| and solving for E yields 



„6+- 



E 



l Ai 



(55) 



Cr~ 



with Ai being an integration constant with respect to 
affinc parameter v. It can be determined by taking the 
limit of (j5"5")) as r tends to some r*. We find 



_ 2n J +2n+l 

Ai = E*r* "~ 2 V '': 



(56) 



We are now able to solve for k using the differential equa- 
tion for E. In fact using (|49p on can write: 



(57) 



K dn 1 I BE 



(1 — k 2 ) dr r E dr 
This equation can be integrated giving 

,r = ! i. - ' 



- 2 e 2 \-2 



(58) 



where A2 is a constant of integration with respect to v. 
To determine A2 we consider that at the point of closest 
approach r — tq (i.e. the closest point that the ray would 



reach in the vicinity of the lensing object), k = 0. As 
consequence from (|58p we have 



A, = 



1 


r ~* +C 




r -2(2r»-l)0>-l) 






3(2n-3) 


m 


r, - 2 


+ Cr.""" 




L ^0 - 2 



-4n* + 10^-7 



(59) 



B. The Deflection Angle 

We now find the form of the deflection angle (f39|) by 
substituing the solutions for E (I5"5"j) and k (|58p giving 



a 



-ct\Jdv — ctQ , 



(60) 



where J is the impact parameter, given in this case by 

i 

3(2ti-3) 



J = v 7 ! - K 2 = 



-4n z + 10n-7 



C 



(61) 



and «i is 



2n z +2n+l 

n-2 



10 11 

i-2 



r* -\- Ct* 



A transformation relation between the affine parameter 
dv and radial distance dr is now needed. We start with 
the propagation equations for <p (|28l) applied to R n grav- 
ity: 

b' = Ek(a-K 

+Ek ( ?^-R + (n - l)R^X (4> + 2A) ) (62) 
in ' 



Substituting in the values of the scalar functions A (|51[) , 
6<\5fy, R ([Ml), X dSU), E §5$ and k ® results in: 
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-2(n + l)(2,i-3) / (2r»-l)(2-. 

»-S r 2 +Cr ^ 



(n - 2)r — — 2(n - 2)r^ - C(4n 2 - 12n + ll)r 



A 2 A? 



(2n 2 - 2n - 1) (4n 2 - lOn + 7) \J + Cr^^ 
Differentiating separately the solution for </> (|52")l with respect to f gives 



, 4 + Sn . „ . in A -In + 11 

2(n - 2)r~ - C (4n 2 - Yin + ll) r — ~ — 



(n-2) 



\ 



(2n 2 -2n-l)(4n 2 -10n+7) 

=—— 7 ~ 4ti 2 -10n + 7 

(n-2) 2 1+CV ^ 



r 



(63) 



(64) 



The transformation relation is obtained by equating 
equations j63) and (HU) giving [H] 



dr = E^a\L 
where L is 



(4.1^-6,1 + 2) 



J 2 r~ 2 + Cr 



(65) 



(n-2) 2 



y (1 + 2n - 2n 2 )(7 - lOn + An 2 ) ' 

Using this transformation relation in (j60[) gives the de- 
flection angle in the form 



(4ti z -C.» - 2; 



J 2 (r- 2 + Cr 



(4n* -1211+11) , 



dr — 7r . (66) 



The standard form of the deflection angle in GR as given 
in [2H | is recovered here when n = 1 in (1661) . 



C. Observables 

We now analyse the behaviour of the deflection angle a 
by computing the deflection angle against n for different 
distances r* from the sources and for different values of 
the distance of closest approach ro. We then plot the 
ratio a to a GR , the GR case corresponding to n = 1, 
against n for the two cases as shown in Fig([rj) and Fig@ 
respectively. As a fiducial system, the distances ro are in 
units of the Schwarzschild radius. 

The divergence of the curves in both plots is indicative 
of the deviation from the standard GR bending angles 
values. In FigJT]), the deflection angle is independent of 
the distance from the source r*. This shows that the 
banding angle does not varies with even if the metric 
is not asymptotically flat. In Fig@ instead one can see 
that for a fixed n, the deflection angle varies for different 
values of tq . In particular for n < 1 there is more bending 



as the values of ro increases whereas for n > 1 more 
bending occurs as ro decreases. 

Fig@ also tells us that, at fixed ro the bending angle 
first decreases with n for n < 1 and then , for n > 1, 
starts increasing. However it is important to remember 
that because of the limits on the validity of the (|5i~]) -([5T 
our conclusions are only valid for n < 1.23. 



VI. CONCLUSION 

In this paper we have used the 1 + 1 + 2 covariant ap- 
proach to derive a general framework for studying strong 
lensing in f(R) gravity. We found that the bending angle 
is dependent of the details of the theory of gravity, but 
this dependence is limited to the feature of the metric 
of the gravitational source. Using the simple example of 
R n gravity and one of its exact spherically symmetric so- 
lutions, we have shown that the banding angle does not 
depend on the position of the observer, in spite of the 
lack of the asymptotic flatness of the metric and that the 
bending angle can be increased of decreased by changing 
the parameters of the theory of gravity. This means that 
using experimental data on the strong lensing one can 
obtain constraints on the function / and consequently 
obtain information on the nature of the gravitational in- 
teraction in the strong field regime. 
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FIG. 1: Plot of the bending angle a compared to the bending 
angle in general relativity against n corresponding to differ- 
ent values of r„, distance from the source. The result shows 
that the bending angle value is independent of r* as the plot 
remains the same when r* is varied. 
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FIG. 2: Plot of the bending angle a compared to the bending 
angle in general relativity against n corresponding to different 
values or ro, distance of closest approach. 
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